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I. INTRODUCTION 



The ultimate goal of high energy heavy-ion studies is to investigate nuclear matter under 
extreme conditions of high density and temperature. In particular, extensive experimental 
and theoretical efforts have been devoted to probe the nuclear Equation of State (EOS), 
to identify the Quantum Chromodynamics (QCD) phase transition and to study properties 
of the Guark-Gluon Plasma (QGP) which is a novel form of matter. New information 
extracted from these studies is critical to many interesting questions in both nuclear physics 
and astrophysics. It is thus fundamentally important to find sensitive experimental probes 
of hot and dense nuclear matter. The collective phenomena found in heavy-ion collisions, 
such as the transverse (sideward directed) and elliptic flow of nuclear matter, have been 
shown to be among the most sensitive probes of the nuclear EOS. The transverse flow is the 
sideward deflection of nuclear matter moving forward and backward in the reaction plane. 
The transverse flow has been observed for nucleons, nuclear fragments and newly produced 
particles, such as pions and kaons, in heavy-ion collisions at all available energies ranging 
from the Fermi energy to the highest RHIC energy Especially, around 4 GeV/nucleon 

where a transition from "squeeze-out" perpendicular to the reaction plane to the in-plane 
flow is expected, a number of flow measurements have been carried out at the AGS/BNL 
I^J^ and the JINR/DUBNA [^-[T^. Clear evidence of the transverse and elliptic flow effects 
for protons and vr" mesons have been obtained in central C-Ne and C-Cu collisions at 
energy of 3.7 GeV/nucleon by SKM-200-GIBS Collaboration at JINR/DUBNA. In C-Ne 
interactions the transverse flow of vr" mesons is in the same direction as for the protons, 
while in C-Cu collisions pions show antiflow behabiour. From the transverse momentum 
and azimuthal distributions of protons and tt~ mesons with respect to the reaction plane 
defined by participant protons, the flow F (a measure of the collective transverse momentum 
transfer in the reaction plane) and the parameter 02 (a measure of the anisotropic emission 
strength) have been extracted. The flow effects increase with the mass of the particle and 
the mass number of target At- The scaled flow Fg = F/{A](^ + A]^^) have been used for 
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comparison of transverse flow results of SKM-200-GIBS with flow data for various energies 
and projectile/target configurations. The Fg demonstrates a common scaling behaviour 
of protons flow values for different energies (Bevalac, GSI/SIS, Dubna, AGS, SPS) and 
systems. Focusing on the total transverse flow, i.e., integrated over transverse momentum, 
these studies have revealed much interesting physics. In the present work, we report results 
of a differential transverse flow analysis as a function of transverse momentum for protons 
and pions in central C-Ne and C-Cu collisions at a beam energy of 3.7 Gev/nucleon measured 
at the SKM-200-GIBS setup of JINR. 

II. EXPERIMENT 

The SKM-200-GIBS setup consists of a 2 m streamer chamber with volume 2 x 1 x 0.6 m^, 
placed in a magnetic field of 0.8 Tesla, and a triggering system. The streamer chamber was 
exposed to a beam of C nuclei accelerated in the synchrophasotron up to a momentum of 
4.5 GeV/c/nucleon (beam energy Ebeam = 3.7 GeV/nucleon). The thickness of the solid 
target Cu (in the form of a thin disc) was 0.2 g/cm^. Neon gas filling of the chamber also 
served as a nuclear target. Photographs of the events were taken using an optical system 
with three objectives. The experimental set-up and the logic of the triggering system are 
presented in Fig. 1. 

The triggering system allowed the selection of "inelastic" and "central" collisions. The 
"inelastic" trigger, consisting of two sets of scintillation counters mounted upstream (Si - S4) 
and downstream (S5, Se) the chamber, has been selecting all inelastic interactions of incident 
nuclei on a target. The "central" triggering system was consisting of the same upstream part 
as in the "inelastic" system and of scintillation veto counters {Sch, Sn), registering a projectile 
and its charged and neutral spectator fragments, in the downstream part. All counters were 
made from the plastic scintillators and worked with photomultipliers PM-30. The Si counter 
with the scintillator of 20x20x0.5 cm^ size, worked in the amplitude regime and identified 
the beam nuclei by the charge. The nuclei from the beam, going to the target, has been 
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selected using the profile counters S2, S3 with the plastic of 15 mm diameter and 3 mm 
thickness and "thin" counter S4 (15 mm and 0.1 mm correspondly). The Sch counters (two 
counters with plastic of 25x25x0.5 cm^ size) were placed at a distance of 6 m downstream 
the target and registered the secondary charged particles, emitted from the target within 
a cone of half angle Gc/i=2.4'^, 2.9*^ (the trigger efficiency was 99% for events with a single 
charged particle in the cone). The S„ counters registered the neutrons, emitted from the 
target in the solid angle 0^=1.8", 2.8*^ (the trigger efficiency was 80% for one neutron). The 
S„ telescope consists from the five counters of 40x40x2 cm^ size, layered by iron blocks 
of 10 cm thickness. The trigger selected central C-Ne and C-Cu collisions defined as those 
without charged projectile spectator fragments (with pz > 3 GeV/c ) within a cone of half 
angle 6^,. = 2.4°. 

The streamer chamber pictures were scanned twice, and a subsequent third scan re- 
solved ambiguities or discrepancies between two scans. Primary results of scanning and 
measurements were biased due to several experimental effects and appropriate corrections 



were introduced [D^OTS 



1) The triggering system selected required projectile nuclei from a primary beam with an 
efficiency higher than 99%. The contamination due to interactions of other projectile nuclei 
with charge less than required biased the values of average multiplicities of secondary par- 
ticles only insignificantly (the correction value < 0.5%). 

2) A trigger bias for central collisions arises whenever a projectile fragment hits the mini- 
mum bias veto-counter system and thus simulates a projectile-nucleus fragment. The effect 
was studied by simulating trajectories of secondary particles generated within the frame- 
work of the cascade model [|1^ . The geometry of the experimental set-up and magnetic field 



distribution were taken into account. The biases thus estimated turned out to be below 
statistical errors in multiplicity, pt and y distributions. 

3) The corrections due to secondary interactions within a solid target turned out to be neg- 
ligible (the correction value < 0.5%). 

4) Pion detection was biased against low momenta. Pions were registered with practically 



no bias when their momenta were J9 > 40 MeV/c for the sohd targets and p > 20 MeV/c for 
gaseous A^e. Pions with lower momenta were also registered, however their detection effi- 
ciency depends on the path length in the target. The corrections could be delivered properly 
only for colliding with equal masses under the requirement of the symmetry y — pt plots in 
the cm. system ( ~ 0.6% in average multiplicity). A very rough estimation for asymmetric 
pairs of colliding nuclei showed that the correction values are (0.5 — 2)%. 

5) The samples of negative secondaries include K~ and S~ tracks. The contamination was 
estimated from pp data and from our results on strange particle production [|l8l (the 
correction value < 0.5%). 

6) Corrections for scanning losses originate from two sources: 

a) scanning inefficiency (~ 2%); 

b) losses of the tracks with a small projection length (which may be screened by the target 
container or a flash around the vertex) and/or a small track curvature ( ~ 1 — 4%) The 
corrections were estimated under the requirement of azimuthal symmetry for each emission 
angle interval for all groups of the measured interactions separately. The correction values 
are 0.5 - 3.5%. 

7) The measurement precision essentially depends on the track length and its dip angle. In 
the samples of events considered ~ 5 — 15% of vr" mesons turned out to be unmeasurable 
or were rejected because of large errors. The measurement losses of tt~ mesons concerned 
practically only well determined intervals of emission angles, and, consequently, appropriate 
corrections, based on azimuthal symmetry, could be introduced in the spectra. Inaccuracy of 
the corrections was calculated, and systematic uncertainties of these corrections were found 
to be less than statistical ones. The ratio o'cent/o'inei (that characterizes the centrality of 
selected events) is (9±1)% for C-Ne and (21±3)% for C-Cu. The average errors in mea- 
suring the momentum and production angle are < Ap/p >= (8-10)% and AO =1° — 2°, 
respectively, for protons, while for pions < Ap/p >= 5% and A© =0.5°, respectively. 
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III. DIFFERENTIAL FLOW OF NUCLEONS AND PIONS IN C-NE AND C-CU 

COLLISIONS 



Several different methods have been utihzed for analyzing flow effects in relativistic 



nuclear collisions, among them the transverse momentum analysis [19| and Fourier expansion 
of particle azimuthal angle distributions pO|j2T[| have been most widely used. More recently, 
it has been found that the differential analysis of the flow strength as a function of transverse 
momentum is more useful in revealing detailed properties of the hot and dense matter, see, 
e.g. P,P^|-pB[]. In this work, we adopt the method of differential flow analysis first used 
by the E877 collaboration in ref. [^]. Let N~^{N~) be the number of particles emitted in 
the same (opposite) direction of the transverse flow near projectile rapidity, then the ratio 
R{pt) = {dN^ /dpt)/{dN~ /dpt) as a function of pt is a direct measure of the strength of 
differential flow near the projectile rapidity. It has been shown experimentally that more 



detailed information about the collective flow can be obtained by studying this ratio p6 



The differential flow data also provides a much more stringent test ground for relativistic 
heavy-ion reaction theories. It was shown theoretically by Li et al in ref. and Voloshin 
in ref. ||2^ that the ratio R{pt) at high pt is particularly useful in studying the EOS of dense 
and hot matter formed in relativistic heavy-ion collisions. 

The data have been analysed event by event using the transverse momentum technique of 
P.Danielewicz and G.Odyniec [^. The reaction plane have been defined for the paricipant 
protons i.e. protons which are not fragments of the projectile {p/z > 3 GeV/c , B < 4°) 
and target {p/z < 0.2 GeV/c). They represent the protons participating in the collision. 
The analysis have been carried out on 723 C-Ne and 663 events. The average multiplicities 
of participant protons are < Np >= 12.4 ± 0.5 in C-Ne and < Np >= 19.5 ± 0.6 in C-Cu, 
respectively. 

The admixture of tt'^ mesons amongst the protons is about (25-^27) %. For the event by 
event analysis it is necessary to perform an identification of 7r+ mesons and separate them 
from the protons. The statistical method have been used for identification of 7r+ mesons. 
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The main assumption is based on the similarity of spectra of vr" and vr"*" mesons (n,r, Pt, 
Pi). The two-dimentional - transverse and longitudinal momentum distribution [pt, pi) have 
been used for identification of mesons. It had been assumed, that 7t~ and vr^ mesons 
hit a given cell of the plane {pt, pi) with equal probability. The difference in multiplicity 
of 7r+ and tt~ in each event was required to be no more than 2. After this procedure the 
admixture of vr^ is less than (5-7)%. The temperature of the identified protons agrees with 



our previous result |29|, obtained by the spectra subtraction. 

In the transverse momentum analysis technique of P.Danielewicz and G.Odyniec [|l9l the 
reaction plane is defined by the transverse vector Q 

n 
i=l 

where i is a particle index, Ui is a weight and p±i is the transverse momentum of particle 
i. The reaction plane is the plane containing Q and the beam axis. The weight uji is taken 
as 1 for ji> Jem and -1 for ji< jcm, where jcm is c.m.s. rapidity and ji is the rapidity 
of particle i. Autocorrelations are removed by calculating Q individually for each particle 
without including that particle into the sum (1). 

^ n 

Qj = ^iP^i (2) 

The transverse momentum of each particle in the estimated reaction plane is calculated as 

E ■ (plj ■ pJa) 

p,/ = ^ (3) 

As we study an asymmetric pair of colliding nuclei, we chose to bypass the difficulties asso- 
ciated with the center-of-mass determination and carried out the analysis in the laboratory 
frame. We have replaced the original weight tUj, by the continuous function uji= Ui- < y > as 
in 



30| , where < ?/ > is the average rapidity, calculated for each event over all the participant 
protons. 

It is known [lin], that the estimated reaction plane differs from the true one, due to the 



finite number of particles in each event. For the estimation of reaction plane resolution we 



divided randomly each event into two equal sub-events (1 and 2) and the reaction plane vec- 
tor in both half-events have been evaluated separately, getting Qi and Q2. Then estimated 
the azimuth angle difference between the vectors Qi and Q2 = $1 — $2- The distribu- 
tions of this have been plotted for the whole experimental samples. The distributions 
are sharply peaked around Ay9^ = ( for C-Cu see Fig. 2). A uniform distributions have 
been recovered for the whole event ($) and two sub-events ($1, $2) reaction planes azimuth 



angles. According to ref. [|T^], the a,. =< Aipji > /2 is a good measure of the reaction plane 
resolution. We found that cr,. = 26° for C-Ne and ar = 23° for C-Cu, respectively. 

First, in Fig. 3 we examine the transverse momentum spectra of protons detected on 
the same {dN~^/dpt) and opposite {dN~ / dpt) side of the transverse flow (reaction plane) in 
central C-Ne (window a) and C-Cu (window b) interactions, respectively. These protons 
are not used in the definition of the reaction plane. Here only protons emitted in the 
rapidity intervals of 1.7< Y < 2.4 for C-Ne and 1.3< Y < 2.7 for C-Cu, respectively, are 
considered. The chosen rapidity ranges are around the projectile rapidity of Yp^oj = 2.28. In 
these rapidity ranges particles with high transverse momenta must have suffered very violent 
collisions and thus originate most likely from the very hot and dense participant region. On 
the other hand, particles with low transverse momenta are mostly from cold spectators. It 
is seen that there is a clear excess of protons emitted to the same side of the directed flow 
for both C-Ne and C-Cu collisions. In both collisions the spectra show typical exponential 
behaviour for p^ > 0.2 GeV/c. The spectra for particles in the same and opposite directions 
of the transverse flow are approximately parallel to each other at pt larger than about 0.7 
GeV/c. These findings are similar to those observed by the E877 collaboration for protons 
in central Au-Au collisions at a beam energy of 10.8 GeV/nucleon p6| . 



The above observations can be understood qualitatively within the transversely moving 



thermal model of Li et al. . Assuming that all or a fraction of particles in a small rapidity 
bin around Yproj are in local thermal equilibrium at a local temperature T, and the center 
of mass of these particles are moving with a velocity (3 along the transverse flow direction 
+x in the reaction plane, then the transverse momentum spectrum of these particles can be 

8 



written as 



C--f{E-P- ptCos{(t))) ■ e-7(^-/3ftco.(<A))/T 



Ptdptd(t)dY 

where C is a normalization constant, 7 = — and is the azimuthal angle with 

respect to the reaction plane. The transverse momentum spectra for particles emitted in 
the same {dN~^ / ptdpt) and opposite {dN~ / ptdpt) directions of the transverse flow are then 

^ = C,.-/-(,ETra)e- (5) 

where a = j^pt/T. These distributions reduce to simple exponents 

^ oc exp{-pt/T^ff) (6) 

at high transverse momenta pt- In the above equation, the inverse slopes or effective tem- 
peratures Tg^j in the semilogarithmic plot of the spectra at high pt are 



—r- — — lim 



d , /rfiV±x 7 
- — - In 



^{cosh{Y)^p). (7) 



.dpt ^Ptdpt' 

In general, the inverse slope ^/T^ff reflects combined effects of the temperature T and the 
transverse flow velocity /3. For the special case of considering particles at midrapidity, (3 
is zero and the effective temperatures are equal to the local temperature T of the thermal 
source. Otherwise, one normally expects T^j^j^ > T~f^. For high energy heavy ion collisions 
in the region of 1 to 10 GeV/nucleon, is much smaller than cosh(Y) around the projectile 
rapidity, thus one again expects T^^j ^ T^jj- at high transverse momenta, i.e two approxi- 
mately parallel spectra dN^ /p^dpt and dN~ /p^dpt at high pt. The distributions {dN^ /dpt) 
and {dN~ /dpt) have been fitted by eq.(5). For C-Ne in the rapidity interval of 1.7< Y < 
2.4 from the {dN^ / dpt) spectra at high Pt { Pt > 0.6 GeV/c) the temperature T and flow 
velocity have been obtained: T = 128 ± 5 MeV/c, /? = 0.016 ± 0.004; from the {dN'/dpt) 
spectra at high pt the following values of temperature T and flow velocity have been obtained: 
T = 130 ± 4 MeV/c, /3 = 0.019 ± 0.005. For C-Cu in the rapidity interval of 1.3< Y < 2.7 
respectively, from the {dN^ /dpt) spectra at highpt ( pt > 0.6 GeV/c): T — 135 ±4 MeV/c, 



(3 = 0.025 ± 0.006; from the {dN-/dpt) spectra: T = 132 ± 6 MeV/c, f3 = 0.023 ± 0.007. 
The results of fitting are superimposed on spectra in Fig. 3. One can see, that the values of 
temperatures and flow velocities from the {dN^/dpt) and {dN~ / dpt) spectra coincide both 
for C-Ne and C-Cu. 

To study the strength of differential transverse flow, we compare in Fig. 4 the ratios 
R{pt) as a function of pt for C-Ne and C-Cu collisions. It is seen that the ratios increase 
gradually at low pt and reach a limiting value of about 1.9 and 2.5 in the reaction of C-Ne 
and C-Cu, respectively. The values of R{pt) are greater than one in the whole transverse 
momentum range, indicating that protons are emitted preferentially in the flow direction 
at all transverse momenta. It is an unambiguous signature of the sideward collective flow. 
Similar results have also been obtained by the E877 collaboration for protons in Au-Au 
collisions at 10.8 GeV/nucleon where the ratio R{pt) increases with pt and flnally saturates 
at about 2. 

The saturation of R{pt) at high pt is also what one expects within the transversely moving 
thermal model. From Eq. ^ one obtains readily 

This ratio normally increases with pt for any given values of T and /?, but it becomes almost a 
constant at high pt for very small /3/T ratios. In relativistic heavy- ion collisions, particles are 
emitted continuously at different freeze-out temperatures during the whole reaction process. 
In particular, particles with high pt are mostly emitted in the early stage of the reaction 
from the most violent regions where the local temperatures are high. Since the ratio R{pt) 
varies very slowly with pt for low 13 /T ratios, one thus expects an approximately constant 
value of R{pt) at high pt. At low transverse momenta, however, R{pt) is affected mostly by 
particles from the cold spectators, and the ratio R{pt) approaches one as pt goes to zero. 

Now we turn to the analysis of differential transverse flow for negative pions in the C-Ne 
and C-Cu reactions. Shown in Fig. 5 are the R{pt) ratios for negative pions in the rapidity 
range of 1.0 < Y < 2.5. The solid lines are the linear flts to the experimental data to 
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guide the eye. In the C-Ne coUisions the ratio R{pt) is greater than one in the whole range 
of transverse momentum, indicating pions are flowing on average to the same direction as 
protons. Moreover, the strength of the pion differential transverse flow increases almost 
linearly with pt and reaches about 1.2 at pt=0.9 GeV/c. On the contrary, it is seen that 
an anti-flow is observed for tt~ mesons in the reaction of C-Cu. Thus, in C-Ne interactions 
the pions are preferentially emitted in the direction of the proton transverse flow, while 
in C-Cu reactions pions are preferentially emitted away from the direction of the proton 
transverse flow. It is necessary to mention here that a similar target dependence of the 
pion preferential emission in asymmetric nucleus-nucleus collisions was first observed by the 



DIOGENE collaboration [^. Our results are consistent with theirs. This target dependence 
of the apparent pion transverse flow was explained quantitatively by nuclear shadowing 
effects of the heavier target within nuclear transport models |]3^j33|l . Our results on the pion 
differential transverse flow and its dependence on the target are thus also understandable. 



IV. SUMMARY 

In summary, the differential transverse flow of protons and pions in central C-Ne and 
C-Cu collisions at a beam energy of 3.7 Gev/nucleon were measured at the SKM-200-GIBS 
setup of JINR. The strength of proton differential transverse flow is found to first increase 
gradually and then saturate with the increasing transverse momentum. /,From the transverse 
momentum distributions of protons emitted in the reaction plane to the same and opposite 
side of the transverse flow the temperatures and flow velocity (3 have been extracted in C- 
Ne and C-Cu collisions. The observations are in qualitative agreement with predictions of 
a transversely moving thermal model. In the whole range of transverse momentum studied, 
pions are found to be preferentially emitted in the same direction of the proton transverse 
flow in the reaction of C-Ne, while an anti-flow of pions is found in the reaction of C- 
Cu due to stronger shadowing effects of the heavier target. The differential flow data for 
both protons and pions provide a more stringent testing ground for relativistic heavy-ion 
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reaction theories. Detailed comparisons with predictions of relativistic transport models on 
the differential flow will be useful to extract more reliable information about the nuclear 
EOS. Such an endeavor is underway and results will be published elsewhere. 
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FIGURE CAPTIONS 



Fig. 1. SKM-200-GIBS experimental set-up. The trigger and the trigger distances are not 
to scale. 

Fig. 2. Reaction plane resolution for C-Cu collisions. The curve is the result of data 
approximation by the polynom of 4-th order to guide the eye. 

Fig. 3. The transverse momentum distributions of protons emitted in the reaction plane 
to the same (o) and opposite (A) side of the transverse flow in a) C-Ne and b) C-Cu 
collisions. The lines are results of fitting by eq.(5) at high pt. 

Fig. 4. Ratios of the yield of protons emitted to the same and opposite side of the transverse 
flow as a function of pt in C-Ne (o) and C-Cu (A) coUisions. The curves are logarithmic 
fits to guide the eye. 

Fig. 5. Ratios of the yield of 7r~ mesons emitted to the same and opposite side of the 
transverse flow as a function of pt in C-Ne (o) and C-Cu (A) collisions. The lines are 
hnear flts to the data to guide the eye. 
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FIGURES 
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FIG. 1. Experimental set-up. The trigger and trigger distances are not scale. 
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FIG. 2. Reaction plane resolution for C-Cu collisions. The curve is the result of data 
imation by the polynom of 4-th order to guide the eye. 
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FIG. 3. The transverse momentum distributions of protons emitted in the reaction plane to 
the same (o) and opposite (A) side of the transverse flow in a) C-Ne and b) C-Cu colUsions. The 
Unes are results of fitting by eq.(6) at high pf. 
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FIG. 4. Ratios of the yield of protons emitted to the same and opposite side of the transverse 
flow as a function of pt in C-Ne (o) and C-Cu (A) coUisions. The curves are logarithmic fits to the 
data to guide the eye. 
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FIG. 5. Ratios of the yield of 7r~ mesons emitted to the same and opposite side of the 
transverse flow as a function of pt in C-Ne (o) and C-Cu (A) collisions. The lines are linear fits to 
the data to guide the eye. 
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